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ABSTRACT 

We calculate the disk and boundary layer luminosities for accreting rapidly rotating 
neutron stars with low magnetic fields in a fully general relativistic manner. Rotation 
increases the disk luminosity and decreases the boundary layer luminosity. A rapid 
rotation of the neutron star substantially modifies these quantities as compared to 
the static limit. For a neutron star rotating close to the centrifugal mass shed hmit, 
the total luminosity has contribution only from the extended disk. For such maxi- 
mal rotation rates, we find that much before the maximum stable gravitational mass 
configuration is reached, there exists a limiting central density, for which particles in 
the innermost stable orbit will be more tightly bound than those at the surface of 
the neutron star. We also calculate the angular velocity profiles of particles in Keple- 
rian orbits around the rapidly rotating neutron star. The results are illustrated for a 
representative set of equation of state models of neutron star matter. 
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1 INTRODUCTION 

For accreting neutron stars in old binary systems, also 
known as low-mass X-ray binaries (LMXBs), a narrow 
boundary layer girdling the neutron star will form next to 
the neutron star surface. The importance of the boundary 
layer derives from the possibility that this could be the site 
for the emission of a variable isothermal blackbody radia- 
tion component observed in the spectra of LMXBs charac- 
terized by very high X-ray luminosity (Mitsuda et al 1984) . 
For weak magnetic field neutron stars, the boundary layer is 
expected to be substantially more X-ray luminous than the 
entire extended accretion disk on general theoretical grounds 
(Sunyaev & Shakura 1986; King 1995). An important fea- 
ture of disk accretion onto a weakly magnetized neutron 
star is that the neutron star will get spun up to its equilib- 
rium period (~ milliseconds), over a timescale of hundreds 
of millions of years (see Bhattacharya & van den Heuvel 
1991 and references therein). A rapid spin of the neutron 
star will enhance its equatorial radius and also relocate the 
inner boundary of the accretion disk closer to the neutron 
star surface. In effect, this would imply a narrowing down of 
the boundary layer separation. Consequently, the boundary 
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layer luminosity is expected to be much smaller in compar- 
ison to the static or slowly rotating neutron star case, and 
this can alter the X-ray emission spectra of LMXBs. 

The efi^ect of rotation of the neutron star on the ac- 
cretion luminosities was considered by Datta, Thampan & 
Wiita (1995), using the 'slow' rotation (but general rela- 
tivistic) formalism due to Hartle & Thorne (1968). These 
authors found that rotation always increases the disk lu- 
minosity, usually decreases the boundary layer luminosity 
and reduces the rate of angular momentum evolution, and 
gave quantitative estimates corresponding to realistic neu- 
tron star models. An important parameter in this connection 
is the radius of the innermost stable circular orbit (rorb)- 
This quantity plays a central role in deciding the magni- 
tude of the gravitational energy release, and hence the accre- 
tion luminosities. The relevance of Varb was emphasized by 
Kluzniak & Wagoner (1985) who pointed out that for non- 
magnetic accreting neutron stars, it is incorrect to make the 
general assumption that the accretion disk will be separated 
from the neutron star surface by a thin boundary layer. The 
boundary layer separation will depend on whether the equa- 
tion of state (EOS) of neutron star matter is stiff or soft. For 
rapidly rotating neutron stars, Cook, Shapiro & Teukolsky 
(1994) calculated the marginally stable circular orbits for 
application to angular momentum evolution of isolated neu- 
tron stars. 
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Accretion onto a rapidly rotating neutron star can bring 
in several interesting features. LMXBs are likely to accrete 
material whose total mass can be a substantial fraction of 
the neutron star mass (^0.1 Mq ). This can severely re- 
duce the magnitude of the boundary layer luminosity (King 
1995). Another important question is whether or not the ac- 
creting neutron star will be disrupted once it reaches equi- 
librium rotation rate with further arrival of the accreted 
plasma. Recently, Bisnovatyi-Kogan (1993) has given a self- 
consistent analytical solution for an accretion disk structure 
around a rapidly rotating non-magnetized neutron star, us- 
ing rigidly rotating polytropic model. This work also gives a 
simple recipe for estimating the accretion luminosities based 
on accreting black hole analogy. In the present paper, we do 
not address the (jucstioii of the disk structure or instability 
at equilibrium rotation rates, but examine how a rapid ro- 
tation rate of the neutron star will affect the boundary layer 
separation and reorder the contribution to the total accre- 
tion luminosity due to the disk and the boundary layer. The 
structure of rotating neutron stars in general relativity are 
calculated using a numerical code developed by us. This code 
is beised on the Komatsu, Eriguchi k, Hachisu (1989) for- 
malism, as modified by Cook, Shapiro & Teukolsky (1994) 
to incorporate realistic neutron star equations of state. This 
formalism is fully general relativistic and is amenable to a 
self-consistent numerical treatment, employing a Newton- 
Raphson type iterative scheme. We find that rotation in- 
creases the disk luminosity and decreases the boundary layer 
luminosity, and for rotation rates near the centrifugal mass 
shed limit, the total luminosity has contribution only from 
the extended disk. Furthermore, for such maximal rotation 
rates, we find that much before the maximum stable gravita- 
tional mass configuration is reached, there exists a limiting 
central density of the neutron star for which particles in the 
innermost stable circular orbit will be more tightly bound 
than those at the surface of the rotating neutron star. We 
also examine the possible modifications in the angular veloc- 
ity profile of the accreted material in Keplerian orbit brought 
on by rapid rotation of the neutron star. 

The format of this paper is arranged as follows. Sec- 
tion 2 gives the formalism and the basic equations to be 
solved. Section 3 outlines the calculation of rotating neu- 
tron star models. The Keplerian angular velocity profiles 
are described in section 4. The results of our calculations 
are summarized in section 5 and a discussion given in sec- 
tion 6. 



2 ACCRETION LUMINOSITIES FOR A 
ROTATING SPACE-TIME 

The space-time around a rotating neutron star can be de- 
scribed in quasi-isotropic coordinates, as a generalization of 
Bardeen's metric (Bardeen 1970): 

ds"" = -e''+^df + e^"{dr^ + r^de^) + 

e'~''r^sine{d(f> - ujdtf (1) 

where the metric potentials 7, p, a, and the angular velocity 
of the stellar fiuid relative to the local inertial frame (u) are 
all functions of the quasi-isotropic radial coordinate (r) and 
the polar angle {6). We use here geometric units: c = 1 = G. 



Since the metric is stationary and ax;isymmetric, the energy 
and angular momentum are constants of motion. Therefore, 
the specific energy E (in units of the rest energy moc^, where 
mo is the rest mass of the accreted particle) and the specific 
angular momentum I (in units of moc) can be identified as 
—po and p3 respectively, where, {/j, = 0,1,2,3), stands 
for the four-momentum of the particle. From the condition 
p^p^ = — 1, we have the equations of motion of the particle 
(confined to the equatorial plane) in this gravitational field 



aT 



(2) 



^2 ^ g2a+7+p _ £,2 f>2 



= E- V^ 



(4) 



Here is the angular velocity of the star as seen by a dis- 
tant observer, dr is the proper time and V is the effective 
potential given by 



V' = 



_7 + P 



1-1- 



e'y-p 



+ 2u3El - up-?. 



(5) 



The conditions for circular orbits, extremum of energy and 
minimum of energy are respectively: 



= 
> 0. 



For marginally stable orbits. 



V,r 



0. 



(6) 

(7) 
(8) 



(9) 



In our notation, a comma followed by one 'r' represents a 
first order partial derivative with respect to r and so on, and 
a tilde over a variable represents the corresponding dimen- 
sionless quantity. 

From the expression for the effective potential and the 
conditions (6), (7) and (9), one obtains three equations in 
three unknowns, namely, r, and I. In principle, if anar 
lytical expressions for e''"'"'', e^™, e^~'' and uj are known, it 
would be a straightforward exercise to solve these equations 
to obtain r, E, and I. In practice, however, this is not so, 
and the solutions for the metric coefficients e^^'' , e^", e^"'', 
and w have to be obtained as arrays of numbers for various 
values of r and 9 using a numerical treatment. Furthermore, 
the condition (9) will introduce second order derivatives of 
7, p, and 00, which means that care has to be exercised in en- 
suring the numerical accuracies of the quantities calculated. 
For this purpose, it is convenient to express E and I in terms 
of the physical velocity v 



{Q, — Lo)re ''sin6 



(10) 



of the stellar matter with respect to a locally nonrotating 
observer (see Bardeen 1972). 

This gives the following expressions: 
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E-ujI = 



p(7+p)/2 



— v 

(7-p)/2 



(11) 



(12) 



Equations (11) and (12) can be recognized as the condition 
for circular orbits. Conditions (7) and (9) yield respectively, 



. ± [e 



+ P,r)+r2(7fr-p?r-)l 



2 + r(7,r - p,r) 



+[2 + r(7,r — p,r)]'y'y,r — e ''ul^r'TV 



2 \ 
■ P,r-) • 







(14) 



where we have made use of Eq. (13) and its derivative with 
respect to r in order to chminate the second order derivatives 
in Eq. (14). The zero of V^rr will give the innermost stable 
circular orbit radius (rorb) and the corresponding v will yield 
E and I. In equation (13), the positive sign refers to the 
co-rotating particles and the negative sign to the counter- 
rotating particles. In this study we have considered only the 
co-rotation case. 

Depending on the EOS and the central density, neutron 
stars can have radii greater than or less than rorb (Datta, 
Thampan & Wiita 1995). The accretion luminosities will, of 
course, be different for these two cases (Kluzniak and Wag- 
oner 1985; Sunyaev and Shatura 1986; Datta, Thampan and 
Wiita 1995). These quantities can be calculated as follows: 
Case (a): Radius of the star (R) greater than rorb. 

If an accretion disk were to form around a relatively 
large neutron star, the ingress of a particle of rest mass mo 
from infinity to the disk boundary (which will be at the 
stellar surface) will release an amount of energy given by: 



Ed 



mo{l - EK{r = R)} 



(15) 



where EK(r — R) is the specific energy of the particle in 
Keplerian orbit at the surface obtained by solving equation 
(13) to obtain vk = v and solving equations (11) and (12) 
with r = R and v = vk to obtain Ik and Exir — R). 

The energy loss in the boundary layer (a very narrow 
gap near the neutron star surface) will be 



Ebl 



mo{EK{r = R)-Eo} 



(16) 



where Eo is the energy of the particle "at rest" on the stellar 
surface (the particle will be moving with the velocity v = v» 

of the stellar fiuid at the surface, where v, is obtained by 
substituting into equation (10) all the relevant parameters 
for r = R) and is calculated by solving equations (11) and 
(12) for E at r = R and v = Vt. 

Case (b): Radius of the star (R) smaller than radius of rorb- 
In this case, the accretion disk will extend inward to a 
radius corresponding to r = rorb- The energy released in the 
disk as the particle comes in from infinity to the innermost 
stable circular orbit will be 



Ed = mo{l-Eorb}- 

The energy released in the boundary layer will be 



Ebl 



mo{Eorb — Eo} 



(17) 



(18) 



1/2 



(13) 



where Eorb is the energy of the particle in innermost stable 
circular orbit, calculated by finding the r — rorb at which 
equation (14) is satisfied and then solving equations (11), 
(12) and (13) for this r to yield Eorb- The energy £^0 of the 
particle on the stellar surface is calculated as described in 



the previous case. 



3 ANGULAR VELOCITY PROFILES 

For slow rotation of the neutron star, the angular velocity of 
the accreted material in Keplerian orbit around it, n(r), will 
have a profile that has a maximum that is located outside 
the neutron star surface. For rapid rotation rates of the star 
(corresponding to angular velocity close to the Keplerian 
value at the surface), a second type of profile for fl{r) is also 
possible, in which S7(r) exhibits no maximum but increases 
monotonically all the way to the surface of the neutron star. 
In such a situation, the accretion torque on the neutron star 
will not be purely advective. It will become possible for the 
viscous torque to transport angular momentum outwards at 
all radii. This can lead to interesting accretion scenarios. 

The Keplerian angular velocity Qk of a particle in an 
orbit around the rotating neutron star is defined as: 



nK{r) 



p(r)v(r) , , 



(19) 



where v is as given in equation (12). The Keplerian angu- 
lar velocity of the particle in an orbit at the surface of the 
neutron star puts a firm upper bound on the angular veloc- 
ity the star can attain (Friedman, Ipser and Parker 1986) 
and hence the boundary layer luminosity when the star at- 
tains this maximum Q. should be zero (Sunyaev and Shakura 
1986). 



4 RAPIDLY ROTATING NEUTRON STAR 
MODELS IN GENERAL RELATIVITY 

Rapidly rotating neutron star structure in general relativ- 
ity for realistic neutron star EOS have been reported by 

Friedman, Ipser & Parker (1986). Their numerical code is 
based on the programmes developed by Butterworth & Ipser 
(1976). Previous models of rapidly rotating neutron stars 
have been based on incompressible fluids and polytropic 
models (Bonazzola & Schneider 1974; Butterworth 1976). 
Komatsu, Eriguchi & Hachisu (1989) have generalized the 
Newtonian self-consistent field method to a general rela- 
tivistic case to obtain structures of rapidly rotating stars, 
again using the polytropic model. This technique was mod- 
ified by Cook, Shapiro & Teukolsky (1994) for realistic neu- 
tron star EOS for purpose of studying quasi-stationary evo- 
lution of isolated neutron stars. A variant of this approach 
based on spectral methods was developed by Bonazzola et 
al. (1993). 
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In this investigation, we have calculated the structure 
of rapidly rotating neutron stars in general relativity using 
a numerical code developed by the present authors, which is 
based on the method due to Komatsu, Eriguchi & Hachisu 
(1989), as modified by Cook, Shapiro & Tcukolsky (1994) so 
as to incorporate realistic neutron star EOS. The results of 
our code agree with the published results of Friedman, Ipser 
& Parker (1986), and results using the code of Stergioulas 
& Friedman (1995) to less than 1%. Also, wherever a com- 
parison was possible, our results agreed with those reported 
in Cook, Shapiro & Teukolsky (1994) to a similar degree 
of accuracy. For non-rotating equilibrium models, we found 
our results to bo within 0.3% of the published results of Ar- 
nett and Bowers (1977). Wc have constructed rmmcrically 
various sequences of neutron stars, starting from the static 
limit all the way upto the rotation rate corresponding to the 
centrifugal mass shed limit. The latter limit corresponds to 
the maximum Q. for which centrifugal forces are able to bal- 
ance the inward gravitational force. Any further increase in 
Cl will lead to disruption of the star. The general relativistic 
expression for this limit can be found in Cook, Shapiro & 
Teukolsky (1994). 

The structure of neutron stars depends sensitively on 
the EOS at high densities. Although the main composition 
of degenerate matter at densities higher than nuclear mat- 
ter density is expected to be dominated by neutrons, sig- 
nificant admixtures of other elementary particles (such as 
pions, kaons and hyperons) are not ruled out. A persistent 
problem in determining the EOS for neutron star interior 
is what to choose for the interaction potential among the 
constituent particles, for which reliable experimental data 
are not available. All calculations involve either extrapola- 
tions from known nuclear matter properties or plausible field 
theoretical approaches using moan field approximation. An- 
other related but unresolved problem is: what is an adequate 
many-body technique to estimate the higher order correla- 
tion terms in expressions for the pressure. In this paper we 
do not address these problems, but choose, for illustrating 
the results of the present study, six EOS models based on 
representative neutron star matter interaction models. This 
is expected to provide a broad realistic set of conclusions. 
A brief description of the EOS models used here is given 
below. 

(A) Pandharipande (1971) (hyperonic matter): One of 
the early attempts to derive nuclear EOS with admixture of 
hyperons is due to Pandharipande (1971), who assumed the 
hyperonic potentials to be similar to the nucloon-nucleon 
potentials, but altered suitably to represent the different 
isospin states. The many-body method adopted is based on 
the variational approach of .Jastrow (1955). The two body 
wave function was taken as satisfying a simplified form of 
the Bethe-Goldstone equation, in which terms representing 
the Pauli exclusion principle were omitted but simulated by 
imposing a 'healing' constraint on the wave function. 

(B) Bethe- Johnson: Bethe & Johnson (1974) devised 
phenomenological potentials for nucleon-nucleon interaction 
that have realistic short-range behaviour. These authors 
then used the lowest order constrained variational method to 
calculate the EOS of neutron star matter. The work of Bethe 
& Johnson (1974) consists of two different parts: (a) deter- 
mination of the EOS for pure neutron gas and (b) derivation 



of a hyperonic equation of state. For purpose of illustration 
here, wc have chosen their EOS model V (neutron matter). 

(C) Walecka (1974) (neutrons): The EOS model of 
Walecka (1974) corresponds to pure neutron matter, and 
is based on a mean-field theory with exchange of scalar and 
(isoscalar) vector mesons representing the nuclear interac- 
tion. 

(D) Wiringa, Fiks & Fabrocini (1988) (UVI4 + UVII): 
These authors gave a model of EOS for dense nuclear and 
neutron matter which includes three-rmcloon interactions. 
This is a non-relativistic approach based on the variational 
method. The three-body potential includes long-range re- 
pulsive parts that arc adjusted to give light nuclei binding 
energies and nuclear matter saturation properties. The au- 
thors have given three models. We consider here their model 
UV14-I-UVII for beta-stable case: neutrons, protons, elec- 
trons and muons. 

(E) Sahu, Basu & Datta (1993): gave a field theoreti- 
cal EOS for neutron-rich matter in beta equilibrium based 
on the chiral sigma model. The model includes an isoscalar 
vector field generated dynamically and reproduces the em- 
pirical values of the nuclear matter saturation density and 
binding energy and also the isospin symmetry coefficient for 
asymmetric nuclear matter. The energy per nucleon of nu- 
clear matter according to Sahu, Basu & Datta (1993) is in 
very good agreement, up to about four times the equilib- 
rium nuclear matter density, with estimates inferred from 
heavy-ion coUison experimental data. 

(F) Baldo, Bombaci & Burgio (1997): have recently 
given a microscopic EOS for asymmetric nuclear matter, de- 
rived from the Brueckner-Bethe-Goldstone many-body the- 
ory with explicit three-body terms. The three-body force 
parameters arc adjusted to give a reasonable saturation 
point for nuclear matter. 

The pressure-density relationship of these above EOS 
models is illustrated in Fig 1. Of the above EOS, model 
(A) is a relatively soft EOS, models (B), (D) and (F) are 
roughly intermediate in stiffness, whereas models (C) and 
(E) are rather stiff EOS. It may be recalled that a stiffer EOS 
will give a higher value for the non-rotating neutron star 
maximum mass and also a higher value for the corresponding 
radius. The composite EOS for the entire span of neutron 
star densities was constructed by joining the selected high 
density EOS to that of Negele & Vautherin (1973) for the 
density range (lO" - 5 x 10^°) g cm~^ , Baym, Pethick & 
Sutherland (1971) for densities down to ~ 10'' g cmr^ and 
Feynman, Metropolis & Teller (1949) for densities less than 
10^ g cm~^ . 

We found that while the results for EOS's (B), (C), 
(E) and (F) were straightforward, EOS models (A) and (D) 
presented the peculiarity of having the maximum S7 model 
as the maximum stable mass model (in agreement with the 
result of Cook, Shapiro & Teukolsky 1994) suggesting that 
EOS's (A) and (D) belong to the Class I EOS and EOS 
models (B), (C), (E) and (F) to Class II as classified in 
Stergioulas & Friedman (1995). 



5 RESULTS 

All the calculated parameters depend on the central density 
(pc) and rotation rate (fi) of the neutron star. In Table 1 we 
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Table 1. The values of the neutron star gravitational mass Mq, disk luminosity Eij, boundary layer luminosity Egj^, the boundary 
layer separation (i.e, the height above stellar surfax;e, where the innermost stable circular orbit is located) h'^ and the Keplerian 
angular velocity Qk for two values of neutron star rotation rates f2 = and Q = Qma and chosen values of central density pc - see 
text for details. The numbers following the letter E in column 9 stand for powers of ten. pc is in units of 10^* g cm~^ ; Hms and 
Clfc are in units of 10* rad s~^. 

EOS Pc Qms M Ed Ebl h+ Qk 

(Mq ) (mo c^) (mo c^) (km) 

f2 = ci = Q,m.s n = Q = cijns n = n = cims n = n = cims n = o n = n, 



(A) 


26.630 


1.105 


1.183 


1.400 


0.057 


0.073 


0.177 


1.578B 




4 


2.052 


0.000 


1.167 


1.106 




*39.121 


1.370 


1.357 


1.585 


0.057 


0.084 


0.246 


1.141B 


- 


4 


4.239 


0.000 


1.018 


1.370 




43.000 


1.434 


1.381 


1.608 


0.057 


0.085 


0.261 


-9.000B 


- 


6 


4.616 


0.223 


1.000 


1.387 


(B) 


10.670 


0.718 


1.133 


1.400 


0.056 


0.058 


0.110 


4.970E 


- 


5 






1.065 


0.718 




-Lo. ( oU 




1 Ann 


1 79*^ 
1. ( 


n np;7 


n n7i 

U.U 1 1 


n 1 Rn 


R 7Sin w 




C 
O 


1 71 Q 


n nnn 
u.uuu 




n S47 




* 22.600 


1.062 


1.694 


2.034 


0.057 


0.084 


0.237 


1.866B 


_ 


4 


5.057 


0.000 


0.815 


1.063 




26.000 


1.119 


1.730 


2.060 


0.057 


0.086 


0.254 


1.160B 


- 


5 


5.619 


0.286 


0.798 


1.082 




35.270 


1.240 


1.757 


2.060 


0.057 


0.088 


0.284 


-9.098B 


- 


4 


6.403 


0.721 


0.786 


1.133 


(C) 


6.363 


0.609 


1.088 


1.400 


0.053 


0.053 


0.090 


3.880B 


- 


5 


— 


— 


0.900 


0.609 




7.413 


0.689 


1.400 


1.818 


0.057 


0.065 


0.128 


1.058B 




4 


0.137 




0.985 


0.689 




*10.610 


0.845 


1.961 


2.493 


0.057 


0.084 


0.217 


2.199B 




4 


5.140 


0.000 


0.704 


0.845 




19.500 


1.042 


2.279 


2.769 


0.057 


0.094 


0.303 


-3.196E 




3 


8.799 


1.408 


0.606 


0.905 




21.750 


1.070 


2.284 


2.757 


0.057 


0.094 


0.311 


-3.877B 




3 


9.020 


1.508 


0.605 


0.918 


(D) 


8.757 


0.700 


1.106 


1.400 


0.056 


0.057 


0.104 


5.593B 




5 






1.033 


0.700 




10.430 


0.803 


1.400 


1.780 


0.057 


0.070 


0.150 


7.670B 




5 


1.297 


0.000 


0.986 


0.803 




*13.530 


0.953 


1.788 


2.239 


0.057 


0.084 


0.223 


1.055E 




4 


4.896 


0.000 


0.772 


0.954 




22.270 


1.217 


2.160 


2.604 


0.057 


0.097 


0.331 


-6.517B 




3 


8.946 


1.640 


0.640 


1.013 


(E) 


3.556 


0.465 


1.059 


1.400 


0.047 


0.047 


0.069 


1.412E 




4 






0.692 


0.466 




4.064 


0.520 


1.400 


1.887 


0.055 


0.057 


0.097 


7.938E 




5 






0.761 


0.520 




* 6.996 


0.674 


2.338 


3.043 


0.057 


0.083 


0.207 


1.836E 




4 


5.641 


0.000 


0.591 


0.674 




11.000 


0.762 


2.572 


3.221 


0.057 


0.089 


0.254 


-4.476B 




4 


8.341 


0.865 


0.537 


0.711 




13.394 


0.794 


2.589 


3.200 


0.057 


0.089 


0.266 


-8.943B 




4 


8.841 


1.038 


0.534 


0.730 


(F) 


9.136 


0.697 


1.106 


1.400 


0.056 


0.057 


0.104 


1.149E 




4 






1.030 


0.697 




11.900 


0.812 


1.398 


1.760 


0.057 


0.070 


0.152 


1.012E 




4 


1.382 


0.000 


0.988 


0.812 




* 20.493 


1.019 


1.732 


2.109 


0.057 


0.084 


0.233 


9.228E 




5 


5.031 


0.000 


0.797 


1.019 




26.000 


1.100 


1.780 


2.135 


0.057 


0.086 


0.257 


-IAS7E 




4 


5.837 


0.380 


0.776 


1.053 




30.890 


1.155 


1.788 


2.126 


0.057 


0.087 


0.270 


-5.163E 




4 


6.199 


0.568 


0.772 


1.080 



have summarized the functional dependences on pc of some 
important parameters in the calculation of luminosities. In 
order to illustrate this dependence, we choose two limits of 
Cl, namely, the non-rotating or static limit {Q, — 0) and the 
centrifugal mass shed limit {Q — Qms), which are the two 
natural ends of a constant density sequence. As expected, 
the functional form of the dependence of structure parame- 
ters on central density for rotation rate at centrifugal mass 
shed limit {fl = flms) is found to be qualitatively similar to 
that at the static limit and so are not explicitly displayed 
here (see, for example Friedman, Ipser & Parker 1986). The 
listed quantities in Table 1 are the values of the neutron 
star gravitational mass Mq, disk luminosity Ed, boundary 
layer luminosity Ebl, the boundary layer separation (i.e, 
the height above stellar surface, where the innermost sta- 
ble circular orbit is located) h'^ and the Keplerian angular 
velocity Qk- Generally speaking, Vorb exhibits three char- 



acteristics: (a) Torf) is non existent, (b) rorb < R, and (c) 
rorb > R- For the first two cases, the Vorb is taken to be the 
Keplerian orbit radius at the surface of the star. The cases 
for which rorb is non-existent are differentiated in Table 1 
by a dashed line under the column for and the cases for 
which Torb ^ R are indicated by the entries for which is 
zero. The central density at which Vorb exactly equals R in 
the rotating case (we call this the radius cross-over density), 
is indicated by an asterisk over the corresponding central 
density. For each EOS model, we choose five central densi- 
ties corresponding to: (i) Mq = 1.4 Mq at f2 = Q.ms, (ii) 
Ma = 1-4 Mq at f2 = 0, (iii) the radius cross-over density, 
at which rorb = Rior Q, = Q.ms, (iv) Mg = maximum mass 
at f2 = Q.ms and (v) Ma ~ maximum mass at f2 = 0. For 
neutron stars corresponding to the EOS models (A) and (D), 
we find that before the maximum stable mass is attained, 
the sequences become unstable to radial perturbations, and 
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Figure 1. Pressure as a function of density for the EOS models (A-F). 



hence the corresponding maximum O model can be taken as 
the maximum mass stable model. Therefore, in Table 1 we 
have listed, for these EOS, parameters corresponding to the 
maximum f2 model instead of the maximum mass model. 

EYom Table 1, it can be seen that typical increases in 
mass (to maintain the central density constant) from the 
static limit is (14-35)%, with the larger changes correspond- 
ing to the lower densities in the stiffer EOS. The correspond- 
ing increase in radius (see Fig. 4) lie in the range (28-44)%. 

Table 1 also shows that the disk luminosity Ed in- 
creases for n = Q,ms with central density whereas, in the 
non-rotating case, it remains constant (except in the low 
central density regime for the stiff EOS for which the Vorb is 
located at the surface of the star). On the other hand, Ebl 
is substantially higher (almost by a factor of 2) than Ed in 
the static limit but becomes almost zero at the centrifugal 
mass shed limit (because of the fact that rotation rate at 
this limit is very nearly equal to the Kepler frequency of a 
particle in the innermost stable circular orbit at the surface). 
Interestingly, Ebl is negative for densities higher than the 
radius cross-over density. 

In Fig. 2, we display the dependence of Q,ma with pc. 
f^ma varies monotonically with pc. Furthermore, seems 
to possess a marked dependence on the EOS; the softer the 
EOS, the larger the value of Q,ma- All curve labels in Figs 2, 
3, 5 and 7 are as indicated in Fig 1. 

The variation of the gravitational mass Ma, baryonic mass 
Mo, stellar radius R, and the radius Vorb of the innermost 
stable circular orbit with respect to for the radius cross 
over density is illustrated in Fig. 3. From Fig 3, it can be 
seen that until O ~ O.Sfims, the structure parameters change 
slowly but for higher rotation rates (f2 ^0.6f2ms), the rate 



of change is more pronounced. These changes are EOS de- 
pendent, being more substantial for stiffer EOS. 

To illustrate how the boundary layer separation varies 
with Pc O and also with the EOS, we give in Fig. 4, plots 
of Torh and R versus pc for two cases of Q,. In Fig. 4, the six 
graphs corresponding to the different EOS models, display 
Rip. = 0) (curve 1), R{n = D,ms (curve 2), rorb(0. = 0) 
(curve 3) and Vorbi^ ~ ^ms) (curve 4). It is immediately 
apparent from the plots that as fl increases from to firns 
for a fixed pc, R (curves 1 and 2) increases and Vort (curves 
3 and 4) decreases. Furthermore, it can be seen that in the 
static case rorb is generally greater than R for the whole 
range of central densities, the exception being for the stiff 
EOS (C) and (E) (in these latter cases, Vorb is less than 
R for a few lower central density values), whereas at the 
centrifugal mass shed limit, rort is greater than R only for 
very high densities. The intersection of curve (2) with curve 
(4) represents the cross-over density pc for which rorb = R 
at the centrifugal mass shed limit. Since at the mass shed 
limit, the stellar rotation rate is almost equal to the angular 
velocity of a particle in Keplerian orbit at the surface of the 
star, the radius cross-over density represents the point at 
which the particle in the innermost stable circular orbit will 
start becoming more bound than a particle at the surface 
of the star. In other words, Ebl will tend to zero at this 
central density and will become negative at some central 
density greater than this. This is indicated by the vertical 
line in the graph. 

In Fig. 5, we plot Ed, Ebl, the total luminosity 
{Ed + Ebl) and the ratio of Ebl to Ed, all as functions 
of Q/Qrns- All these plots are for radius cross-over density 
for the various EOS. Among the various EOS, Ed differs 
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Figure 2. The neutron star rotation rate at the centrifugal mass shed limit as a function of central density for various EOS models. The 
labels in Figs 2, 3, 5 & 7 are as in Fig 1. 
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Figure 3. Neutron star gravitational mass (Mq), baryonic mass (Mq), radius (R), and the radius {rorb) of the innermost stable circular 
orbit as a function of the star's rotation rate for various EOS. The graphs correspond to the radius cross over central density for each 
EOS. 
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Figure 4. Neutron star radius (i?) and radius (rarb) of the innermost stable circular orbit as functions of central density for two values 
of stellar rotation rate (fl). See text for the details. 



most when CI is large. In contrast, Ebl differs (among the 
various EOS), most for slow rotation rates. Therefore, the 
total luminosity follows a variation similar to that of Ebl- 

The dependence of Qk on Q is shown in Fig 6. For 
convenience of display and comparison, these are taken in 
units of ClmB, and two chosen values of pc- In Fig 6 curve 1 
refers to pc corresponding to Mq = 1.4 Mq at static limit 
and curve 2 refers to Pc for the radius cross-over density. The 
qualitative differences in the behaviour of curve 1 for all the 
graphs is quite noticeable. These differences arise entirely 
due to the differences in the location of the innermost stable 
circular orbit with repect to the stellar surface. 

In Fig. 7 we give plots of the angular velocity profiles 
for the radius cross over densities for each of the EOS. 



assumption made is that the magnetic field of the neutron 
star is too small to affect the accretion flow. It is relevant 
to ask if a quantitative estimate is possible of how low the 
magnetic field should be for the validity of our calculations. 
The Alfven radius (va) is defined by the relationship (see 
Lamb, Pethick & Pines 1973) 

= pivAyivA) (20) 

where p and v are respectively the density and radial ve- 
locity in the accretion disk. The Alfven radius determines 
the location at which magnetic pressure channels the flow 
from a disk into an accretion column structure above the 
magnetic poles. Lamb, Pethick & Pines (1973) show that 



6 DISCUSSIONS 

In this paper we have investigated in a general relativistic 
manner, the effect of rapid rotation on the boundary layer 
and disk luminosity for an accreting, old neutron star. The 



VA 



2.6 X 10'' 



M^/^(M/Mo )V^ ' 

^37 -"-6 



(21) 



where /ttso = BqR^ /lO^" G cm^ , L37 is the total luminosity 
in units of 10^'^ ergs , Re = R/W^ cm and So is the 
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Figure 5. The accretion luminosities in units of the rest mass of the accreted particle as functions of neutron star rotation rate. The 
graphs correspond to the radius cross-over central density for each EOS. 



magnetic field on the surface of the neutron star in gauss. 
The condition that va < R implies that (for the reasonable 
choice: M = 1.4 Mq and Re = 1): 

So < 5.5 X IO^lK^ (22) 

and is necessary for the scenario we have considered to be 
fully self-consistent. In our notation, L = {Ed + Ebl)M(? ^ 
with M the mass accretion rate. 

The main conclusions of this paper are that rotational 
effects of general relativity increase the disk luminosity and 
more importantly from the point of observations, decrease 
the boundary layer luminosity. These effects are small in 
magnitude for small values of Q, but become substantial for 
rapid rotation rates of the neutron star. The boundary layer 
luminosity becomes inconsequential for rotation rates near 
the centrifugal mass shed limit. For such cases, the role of 
radiation pressure on the accretion flow (see Miller & Lamb 
1996) must be re-examined. Also for such cases, accretion 
induced changes in the surface properties of the neutron star 
is an important question to investigate. The vanishing of the 
boundary layer luminosity for rapid rotation rates as found 
in this study is not apparent in a similar calculation using 
the "slow" rotation approximation based on the Hartle & 
Thorne (1968) metric (see Datta, Thampan & Wiita 1995). 
The total luminosity remains fairly constant upto rotation 
rate of about Q.QQ,ms, but declines rapidly to the value of 
the disk luminosity for higher rotation rates. We have not 
considered here the angular momentum evolution of the ac- 
creting neutron star, but have calculated the accretion lumi- 
nosities for chosen fixed values of Q,. So, the boundary layer 



luminosity values listed in Table 1 do not include correc- 
tions for the energy that may go into the spinning up of the 
neutron star. We have not considered the effect of viscosity 

on the accretion luminosity. The viscosity effects require a 
full hydrodynainic treatment and also the radiative transfer 
phenomenon. 

An interesting conclusion of the present study (see Fig 
4) is that neutron star configurations with high central den- 
sities have their innermost stable circular orbit located exte- 
rior to the star. For such configurations that are rotating at 
the centrifugal mass shed limit, particles in the innermost 
stable circular orbit are more bound than particles at the 
surface of the star. This could lead to the formation of an 
inner disk torus. The idea of an inner disk torus has been 
invoked as a possible explanation of flaring branch phenom- 
ena observed in certain Quasi Periodic Oscillators (Kuulkers 
& van der Klis 1995), with radiation pressure playing an key 
dynamical role. Our study seems to suggest that an inner 
disk torus can be formed in the absence of a substantial radi- 
ation pressure, purely as a consequence of general relativistic 
rotational space-time in situations where the rotation rate 
of the accreting neutron star is close to the centrifugal mass 
shed limit. 
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Figure 6. The Keplcrian angular velocity of a particle in the innermost stable circular orbit as a function of the neutron star rotation 
rate. Curve 1 refers to the central density corresponding to a gravitational mass of 1.4 Mq in the static limit, and curve 2 corresponds 
to that at the radius cross-over central density. 
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